Generic Model of Morphological Changes in Growing Colonies of Fungi 
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Fungal colonies are able to exhibit different morphologies depending on the enviromental con- 
ditions. This allows them to cope with and adapt to external changes. When grown in solid or 
semi-solid media the bulk of the colony is compact and several morphological transitions have been 
reported to occur as the external conditions are varied. Here we show how a unified simple mathe- 
matical model, which includes the effect of the accumulation of toxic metabolites, can account for the 
morphological changes observed. Our numerical results are in excellent agreement with experiments 
carried out with the fungus Aspergillus oryzae on solid agar. 



I 

C 

O 
o 



> 

O 

in 
o 



o 

c3 



-a 
c 

o 
o 



X 



I. INTRODUCTION 

Since the pioneering work of D'Arcy Thompson JjJ at 
the beginning of this century understanding the origin of 
complex structures involving living organisms has been 
a challenge for physicists and biologists alike. In par- 
ticular the patterns formed in growing colonies of mi- 
croorganisms, like bacteria J|-0| an d fungi f^-|l6f| have 
attracted much interest because they exhibit properties 
akin to those in inanimate systems far from equilibrium 
p7[ . Early mathematical models of fungal growth ]l8],|l9| , 
considered fungi as purely additive assemblages of dis- 
crete individual hyphal units that duplicate at regular 
intervals and homogeneously. These models predict the 
crossover from exponential growth to linear growth of the 
peripheral zone as the colony matures. However, they 
are unable to account for many of the observable macro- 
scopic properties (as morphological changes for instance) 
exhibited by fungi. 

Fungal growth experiments on agar plates have shown 
[p| |l6| that filamental growth leads to a fractal colony 
with patterns that have much in common with those ob- 
served in diffusion-limited aggregation (DLA) processes 
Similar branched morphologies have also been 



where the relevance of 
|j. Sophisticated corn- 



observed in bacterial growth |!H5 
DLA was nicely demonstrated 

munication mechanisms, including chemotactic signals 
and feedback response, were later proposed || in order 
to understand the origin of branched growth in bacteria. 
A similar mathematical modelling might be adequate to 
explain the DLA- like patterns of hyphae. However, quan- 
titative measurements of the fractal dimension of mycelia 
p3|,p[ have shown that the dimension varies in time and 
from one species to another. Also, the growth mechanism 
of DLA clusters seems rather different from the hyphal 
extension of fungi. Some reaction-diffusion models have 
also been studied |^2|,^3| , which seem to capture esential 
features of the filamental fungi. 

In contrast, under appropriate environmental condi- 
tions, colonies of many fungal species are also able to 



display a compact phase [£4|,^5|] with uni-cellular non- 
cooperative growth and generate non-DLA patterns. In 
this case, transitions between different morphologies are 
also observed [[l4]-[l6) when the relevant parameters in 
the growth conditions are varied. A mathematical model 
able to include all these fungal patterns on solid media 
has remained elusive up to now |l6|]26| ]. Here we show 
that a simple stochastic model containing some generic 
biological mechanisms is able to reproduce a variety of 
patterns observed in the growth of fungal colonies on 
solid agar plates. The model takes into account the ef- 
fect of the accumulation of waste products, which play 
an important role as inhibitors for cell division. 

Matsuura and Miyazima Jl4||l5| carried out experi- 
ments with the fungus Aspergillus oryzae on agar plates 
of different stiffness and with different concentrations 
of nutrients (glucose). They observed that the fun- 
gal colonies on solid media formed a variety of shapes 
(see Figure 1). In nutrient-rich conditions, the colonics 
showed a homogeneous and smooth growth front (Fig. 
la). In contrast, growth in glucose-poor conditions re- 
sulted in colonies with rougher surfaces (Fig. lb). The 
physiological ability of the fungus to divide and grow 
is diminished at lower temperatures |24| , p5[ ], when even 
morphological instabilities (Fig. lc) or an incipient 
branched growth (Fig. Id) can be observed depending 
on the nutrient concentration. Similar patterns have re- 
cently been obtained in experiments with colonies of the 
yeast Pichia membranaefaciens on solidified agarose |l6| . 
How can these morphological changes be qualitatively de- 
scribed and understood within one simple mathematical 
model?. 



II. THE MODEL 

It is known that metabolic products generated by 
fungal cells during the degradation of nutrients diffuse 
through the medium and may act as impurities inhibit- 
ing further cell division 24 2^]ll|. We have considered 
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the effect of toxic metabolites on the morphology of the 
colony. We have studied a model in which the growth 
probability at every site is coupled with an inhibitory 
field generated by the fungus itself. The model is defined 
on a two dimensional lattice as follows. At time t a lat- 
tice site can be either occupied by a fungal cell or vacant. 
Two quantities are assigned to every site x in the lattice: 
(1) the total age of the nearest neighbour cells A(x, t) 
of that site and (2) the field c(S, t) containing the con- 
centration of waste products. The spread of inhibitors 
from a single cell situated at x is modelled by means of 
a diffusion equation with a source: 

= D V 2 c(f , t) +s6(x- So) Q(t - r ), (1) 

where is the step function (0(u) = 1 for u > and 
Q(u) =0 for u < 0). D is the diffusion constant of the 
metabolites, s is the inhibitor chemicals' production rate, 
and To is the time at which the polluting cell at So was 
born. The total density of waste products at any vacant 
site at time t is calculated by adding the contributions 
coming from every occupied site in colony situated within 
a distance d. 

Colony growth occurs because of the division of in- 
dividual cells, thus only nearest neighbours of occupied 
sites have a chance of becoming occupied. We assume 
that the probability for a vacant site S of being occu- 
pied in the next time step increases with the total age 
A(x, t) of the occupied nearest neighbours of that vacant 
site x. Therefore, at every time step, we assign a growth 
probability P(x, t) to every vacant site S given by 

P(S,t) = F[9 A(S,t) <j>(c)], (2) 

which depends on the two fields A(x, t) (total age of 
neighbouring occupied sites) and c(x, t) (waste concen- 
tration) introduced above. In Eq.(Q), F is an arbitrary 
monotonous increasing function satisfying F(0) = and 
-F(oo) = 1 (properties required for P to be a probability). 
The results do not depend on the detailed form of F and 
in the following we use F(x) — tanh(x) for simplicity. 
The suppression of cell division due to the accumulation 
of metabolites is represented by 4>{c) that is a decreasing 
function of the concentration of waste chemicals. In or- 
der to see a significant effect of the inhibitors, the decay 
of <j) with c must be rapid enough. In the simulations we 
are presenting in this paper <j)(c) — exp(— c/co) was as- 
sumed, Co being a concentration threshold above which 
the effect of inhibitors becomes important. Other sim- 
ple functional forms may be used and similar patterns 
are obtained provided that 0(c) decays rapidly enough. 
Finally, the external parameter 9 controls directly the 
growth rate and is associated with the food supply (me- 
dia richer in nutrients correspond to larger values of 9). 



III. RESULTS 

We carried out extensive simulations of the model. The 
simulations were performed on a two dimensional trian- 
gular lattice with open boundary conditions, simultane- 
ously updating Eqs.(0) and (||) in every time step. The 
position of the front is determined and the height h(x, t) 
at every column x is used to study the front dynam- 
ics. In order to mimic the experiments, the growth is 
initiated from a straight line of randomly occupied sites 
at h(x, 0) = 0. The dynamics of the growing front is 
characterized by the height fluctuations, which are mea- 
sured by means of the front width over the total system 
of size L, W(L, t) = y/((h - (h)) 2 ). The use of other 
quantities to measure the width is helpful when instabil- 
ities or overhangs appear, like for instance the difference 
Max[h] — Min[h] between the maximum and minimum 
heights, and gave similar results. 

The model has two external parameters, namely: the 
growth rate 9, and the concentration threshold Co for in- 
hibitory effects. These are to be thought of as effective 
parameters and their relation to experimental variables 
such as temperature, concentration of nutrients, etc. may 
not be a simple one. It is to be noted that, if experimental 
conditions change, for instance in such a way that the ab- 
sorption of inhibitors by the fungal cells is enhanced, the 
effect may be well represented by either a lower threshold 
Co or a wider spread distance d (or even a combination 
of both). 

We performed computer simulations of the model for 
different values of 9 and Co. In the following we discuss 
the different front morphologies that are obtained de- 
pending on the values of these two parameters, and how 
they are related to experimental observations. In Fig. 2 
we present a summary of the resulting numerical fungal 
patterns. 



A. High threshold for inhibitory effects 

The first morphological transition occurs in the limit 
Co — * oo, in which the inhibitors become irrelevant 
(4>(c) w 1). We found that the surface of the colony 
is flat (Fig. 2a) when inhibitors have little effect (for 
Co ^> 1) and the medium is rich in nutrients. Conversely, 
the growth front becomes rougher and Eden-like (Fig. 
2b) for a poorer medium. In fact, a quantitative measure 
of the roughness (see below) indicates that the fractal 
dimension of the front is 1.5 |g6[ as in the Eden model 
p7[ . These morphological phases are separated by a con- 
tinuous phase transition at a critical value of 9 = 9 C . 
For large growth rates (rich medium) 9 > 9 C the front 
is flat, W ~ const. In contrast, if the growth rate is re- 
duced below 9 C (poor medium) the front becomes rougher 
and exhibits the dynamical scaling behaviour typical of 
a scale- free roughening process |2^,^J ■ As can be seen in 
Fig. 3, in the rough phase the width grows like a power 
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law W ~ t 13 until a stationary regime is reached and 
W ~ £ Q , where L is the system size. Our determination 
of the time exponent (3 = 0.24 ± 0.02 and the roughness 
exponent a = 0.46 ± 0.05 indicates that the model be- 
longs to the Edwards- Wilkinson universality class J2l| in 
this phase. The computation of the roughness exponent 
a allows us to determine the fractal dimension Dj of the 
front through the relation Df — 2 — a ~ 1.54, which 
is consistent with Eden-like growth. This morphological 
transition from rough to flat growth has recently been 
studied in detail by us pq | for a much simpler version of 
this model, which corresponds to the limit Co — > oo. We 
found that 9 C = 0.183 ±0.003 in the limit case of Co = oo. 
When approaching the roughening transition from the 
rough phase (9 < Q c ) to the flat phase (6 > 9 C ), a diverg- 
ing correlation lenght £ ~ \9 — 6 c \~ L ' appears. The critical 
exponent was found to be v = 1.10. At the phase transi- 
tion the scaling exponents can be calculated by m app ing 
the problem onto directed percolation (see Ref. (26| for 
further details). Experimental realizations of transitions 
from a flat to a rough front can been observed in the Mat- 
suura and Miyazima experiments shown in Figure la-lb. 
As our model predicts the richer the medium the Hater 
the front. 



B. Low threshold for inhibitory effects 

If the effect of toxic inhibitors on the fungus is en- 
hanced, further interfacial instabilities in the front ap- 
pear. Lower temperatures are less favorable for the phys- 
iological activity of the fungus |24|,|2^] and its capacity to 
divide. One can model this effect by a smaller threshold 
Co for the inhibition of cell division. One can observe 
(Fig. 2c) that, as Co is decreased, deep grooves and over- 
hangs develop in the colony and some regions with high 
concentration of inhibitors are eventually surrounded by 
the colony front. The presence of a unstable growth dy- 
namics can be realized by the appearance of an extremely 
rapid (exponential) growth of the interface width W(L, t) 
in this regime, as shown in Fig. 4. This regime is to be 
compared with the patterns observed in the expriments 
(Fig. lc). 

Finally, in the worst of the situations, the growth be- 
comes very localized and only the dominant filaments 
branch out into the empty medium (Fig. 2d). This is 
due to the building up of extremely high concentrations 
of waste products in some regions that the fungus cannot 
invade any longer. There is qualitative agreement with 
the experiment shown in Fig. Id. 

The four major morphological phases are summarized 
in the qualitative phase diagram in Fig. 5. The solid line 
that separates the rough from the flat phase in Fig. 5 has 
recently been shown to be a nonequilibrium continuous 
phase transition and can be mapped into directed per- 
colation pq |. However, we were unable to find a sharp 
transition between the remaining phases. 



IV. DISCUSSION 



Our results indicate that the morphologies generated 
during the growth of fungal colonies on solid agar can 
be modeled by taking into account (1) the building up 
of a concentration of waste inhibitors in the medium and 
(2) the provision of nutrients. We have studied a simple 
model which allows us to understand the interplay be- 
tween these two mechanisms. The morphological tranfor- 
mations that we have found can be compared with recent 
experiments on pattern formation in fungal growth. 

In our model, a constant growth rate 9 has been as- 
sumed. The effect of the exhaustion and diffusion of nu- 
trients can be modeled by replacing 9 by a diffusive field, 
which is depleted at a certain rate at sites occupied by 
fungal cells. This would have a major impact on the 
front dynamics as the front will eventually stop growing 
because of agar exhaustion. 

Our model approach to pattern formation in fungal 
colony growth is stochastic and includes in a natural way 
the random nature of cell division. This together with 
the accumulation of inhibiting chemicals produces some 
of the different morphological phases observed in experi- 
ments. This is in contrast to other existing models based 
upon non-linear diffusion, in which evolution is determin- 
istic and patterns arise from nonlinear instabilities. They 
represent a different level of description of the same phe- 
nomena. Both type of models coexist at different levels 
of description. 

Finally, while our model can account for the morpho- 
logical changes observed in some fungi growing in agar 
culture, we should mention that this model may not 
be appropiate for describing the morphological patterns 
observed in other fungal species. Our results compare 
well with the fungus Aspergillus oryzae which, like many 
members of this genus (and also the genus Penicillum) , 
form dense colonies that expand only very slowly and 
seem to be self-inhibitory. 
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viding us photographs of their experiments. We thank T. 
Sams for useful discussions on his work and M.J. Carlile, 
S. Matsuura, and E. Calle for helpful comments and cor- 
respondence. This work has been supported by a grant 
from the European Commission (Contract No. HPMF- 
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FIG. 1. Observed patterns of colonies of Aspergillus oryzae 
grown on solid medium (1.5 wt% of a Czapek synthetic agar 
in 25 ml of sterile medium). Experiments were carried out 
in media with two different concentrations of nutrients: 0.1 
wt% glucose for a nutrient-rich medium and 0.01 wt% glucose 
in nutrient-poor conditions, (a) Smooth colony 10 days after 
inoculation at 24°C in a nutrient-rich medium; (b) patterns 
8 days after inoculation at 24° C in a poor medium, where 
the colony developed a rough front, (c) 30 days at 18°C in 
a rich medium. Note the tip splitting dynamics and groovy 
structure of the colony associated with the existence of a mor- 
phological instability, (d) Prominent filamental patterns of 
a colony grown in a nutrient-poor medium after 35 days at 
18°C. Photographs kindly provided by S. Matsuura. 
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FIG. 2. Fungal patterns obtained by numerical simulation 
of the model on a triangular lattice of lateral size N = 500. 
Each front corresponds to the position of the colony sur- 
face at equal time intervals. When the threshold Co is large 
(co = 5.0), two different morphologies appear: (a) the colony 
is flat in a nutrients-rich medium, 6 = 0.1 and (b) rough in 
the nutrients-poor case, 6 = 0.01. For small values of the 
threshold Co, the colony exhibits a morphological instability 
as in (c) for Co = 1.0 ; or prominent filamental growth as 
in (d) for c = 0.7. In the simulations shown, the produc- 
tion rate of toxic metabolites and the diffusion constant are 
fixed to s = 0.05 and D = 1 respectively. Sites within a dis- 
tance d = 16 are considered to calculate the concentration of 
inhibitors at any site. 
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FIG. 3. Interface width of the front vs. time in 
the rough phase (9 — 0.01) for different system sizes 
L = 16, 25, 32, 50, 64, 75, 100 . The slope of the dotted line 
corresponds to the time exponent /3 = 0.24±0.02. In the inset 
the values of the width at saturation are plotted vs. system 
size. The dotted line fits the data and gives the roughness 
exponent a = 0.46 ± 0.05. 
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FIG. 4. Time behaviour of the interface width in the rough 
(solid line) and unstable phases (dashed line). Note the 
change from power law bahaviour W(t) ~ t 13 to exponential 
growth W(t) ~ exp t in the unstable case. The lines are to 
guide the eye and have slopes 0.25 and 1. The system size was 
L = 500. The external parameters are 6 — 0.01 and cq = 5 
(rough phase) and 9 = 0.01 and c = 1 (grooves' phase), 
corresponding to Figs. 2b and 2c respectively. 
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FIG. 5. Sketch of the phase diagram showing the different 
morphologies observed in the model. The solid line corre- 
sponds to the continuous phase transition studied in Ref. [26] . 
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